Multiple Electromagnetic Excitation in Fast 
Peripheral Heavy Ions Collisions 

S. Typel * 

Physique Nucleaire Theorique et Physique Mathematique, 
Universite Libre de Bruxelles, Campus de la Plaine CP229, 
Bvd. du Triomphe, B-1050 Bruxehes, Belgium 

G. Baur 

Institut fiir Kernphysik, Forschungszentrum Jiilich GmbH, 
D-52425 Jiilich, Germany 



*E-Mail: stypel@ulb.ac.be 



1 



Abstract 

We study the corrections of first order electromagnetic excitation due 
to higher order electromagnetic interactions. An effective operator is intro- 
duced which takes these effects into account in the sudden approximation. 
Evaluating the matrix-elements of this operator between the relevant states 
corrections to the first order result are obtained in a simple way. As an 
example we discuss the excitation of the first excited state in ^^Be. It tends 
to improve the agreement between experiment and theory. 
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Electromagnetic Excitation in the energy domain of several tens of MeV/u up 
to relativistic energies is a growing field of study. The cross-section can become 
large and irreducible nuclear effects can be kept under control. With increasing 
beam energy the equivalent photon spectrum becomes harder, and also particle- 
unstable states can be reached. The Coulomb dissociation ^^Li ^Li + 2n and 
^^O — * ^^N + p, which is also astrophysically relevant, are examples Recently, 
bound states were also excited and their (Doppler shifted) de-excitation 7-rays 
were measured. A large deformation of the neutron-rich nucleus '^^Mg was re- 
cently deduced from a measurement of the 2+ —>■ 0^ 885 keV transition to the 
ground state after medium energy electromagnetic excitation [^. The 320 keV 
I -^1^ 7-transition in ^^Be was recently observed. The measured cross-section 
for the ^^Be | ) Coulomb excitation was found to be noticeably less than 

expected from the known lifetime and 1^* order pure Coulomb excitation ||^. 
Apart from possible nuclear and Coulomb-nuclear interference effects, a possible 
reason for this discrepancy is the influence of higher order electromagnetic inter- 
action. It is the purpose of this letter to describe a framework suitable for fast 
projectiles. E.g. the rather loosely bound ^^Be in its 1^* excited | state could 
easily be excited electromagnetically into the continuum in a second step |Q. 

Electromagnetic excitation is mainly characterized by two parameters, the 
adiabaticity parameter 



ujb 

'JV 
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and the strength parameter 



jEX) _ Ze{f\M{EXf,m 



The excitation energy is given by huj, the impact parameter in a straight-hne 



approximation is denoted by 6, and 7 = y 1 — (f/c)^, where v is the projectile 
velocity. The target charge number is denoted by Z, and J^{E\ii) the electric 
multipole operator. 

In Coulomb excitation below the barrier, multiple electromagnetic excitation 
is usually treated in a coupled channels approach using the relevant states from 
appropriate nuclear models, like the harmonic vibrator or rigid rotor. For a review 
see Ref. 0. The situation for Coulomb excitation above the barrier becomes 
simpler, because the excitations tend to be sudden. While Tiu is restricted to a 
few MeV, the adiabaticity parameter ^ is typically less than 1, for the important 
range of impact parameters b > Ri + R2, where Ri and R2 are the nuclear radii 
of projectile and target. Thus fast collisions become the domain of the sudden 
approximation , or of a recently developed low-^ approximation @, |^ . In this 
case it can be advantageous to construct operators which take into account the 
influence of intermediate states. 

For simplicity, let us use the straight-line approximation and the dipole ap- 
proximation. The first order excitation amplitude is given by 

afhEl,0 = -^m■{m^) (3) 



with 



hvb 
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V 
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where the projectile moves in the z-direction and the impact parameter points to 



(5) 



the x-direction. The dipole effective charge is given by 

(1) _ Zbiric - Zjrib 

using a model of two pointlike inert clusters b and c with charge numbers Zf,, Zc 
and masses nib and nic- Kq and Ki are the modified Bessel functions. For ^ — > 
we get the classical Coulomb push 



(6) 



In the limit ^ << 1, the excitation amplitude can be evaluated easily to all orders 



in the sudden approximation. It is given by 



a}f'*""(El) = {f\e 



(7) 



This could be generalized to higher multipolarities and trajectories corrected for 
Coulomb deflection. By comparing a^j}{El^^ — 0) and afl;'^'^^"'(El) the influence 
of multiple electromagnetic excitation can be assessed. This is a remarkably 
simple procedure, all intermediate states are included. The reduction of the 



excitation probability due to higher order effects is given by 

„suaaeni rp-w ^ 

r 
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(8) 



For loosely bound states, e.g. in ^^Be = ^'^Be + n, we use simple model wave 
functions to reveal the characteristic parameters. We choose two models, with the 
correct asymptotic behaviour. Their differences give a feeling about the model 
dependence. For the wave functions we make the ansatz for the radial part in 
the initial state (with orbital angular momentum li = 0) 

/ X / exp(— a,r) , , 

r 

which corresponds to the solution of the Schrodinger equation with a 5-like po- 
tential. For the final state (// = !) we choose 

,4(r) = v^5^^eh^ (10) 

and the more extended wave function 



Uf {r) = 2yaj exp(— a/r), (11) 

respectively. The constants {z = i, f) are calculated from the binding energies 

= With these wave functions the calculated mean lifetimes of the k 

state are 110 fs and 97 fs, respectively. They are smaller than the experimental 
value of (166 ± 15) fs P]. Despite this difference, the ratio of higher order to 
first order effects can be given with some confidence, where, e.g., common spec- 
troscopic factors cancel out. The absolute value depends on more sophisticated 
details of the nuclear model (see, e.g., Ref. 0). The reduction factor r can be 
calculated analytically. We get 

= (^^{l- -arctan^;^^ (12) 



3 / z 



and 

resp., where we have introduced the parameter 

(14) 

The parameter z is directly related to X/f^"* (Eq- 0)- It describes the ratio of the 

strength of the Coulomb push q and the "looseliness" of the system, and is a 

measure of the importance of higher order effects. The quantity r is plotted in 

Fig. 1 for the two model wave-functions described above. 

For large b the sudden approximation fails oo), on the other hand, higher 

( F^ ^ 

order effects diminish due to the decrease of the strength parameter Xfi ■ The 
product 

^Ei) , Ze{f\M{Elf^)\t)uj 

is a very small number for low excitation energies Tluj and high projectile velocities 
V. The ranges of validity for the first order calculation {x small, ^ arbitrary) 
and the sudden approximation {S, small, x arbitrary) overlap. In a convenient 
and accurate interpolation procedure we calculate the total cross section in the 
following way 



a(~) = 27r/ r(h)\a^}}{El,i(h))\'^hdh (16) 



which should be an accurate expression for all values of b in the integrand. For 
small impact parameters we have ^ ~ and the first order approximation cancels 



out in the calculation of a^°°\ We compare it with the total cross section in the 
first order calculation 



Both cross sections depend on the value of the minimum impact parameter bmin- 
The change in the two cross sections will be similar so that their ratio is less 
affected by a change in bmin- If the sudden approximation is not well enough 
fulfilled, one could use the low-,^ approximation ||^, ^ which takes second order 
electromagnetic effects into account. 

For the ^^Be Coulomb excitation a reduction of the cross section from (490 ± 
50) mb (expected from the first order calculation) to the measured (191 ±26) mb 
was recently found in an experiment at GANIL Indeed, the parameter z 
can become substantial in this case, and higher order effects are not negligible. 
For a collision with 6 = 15 fm and an energy of 45 ■ A MeV (f/c ~ 0.32) we 
have z ~ 0.38 and from Fig. 1 we see a substantial reduction of the excitation 
probability. In a recent experiment at RIKEN 0, which is currently evaluated, 
an energy of about 65 ■ A MeV {v/c ~ 0.38) for the ^^Be was used. This leads 
to a value of z ~ 0.32 for the same impact parameter corresponding to a smaller 
reduction. 

The product X/f^^ ' ^ takes on the small values 0.00375( JjMjl/i| JjMj) and 
0.00399( JiMjl/i| J/M/), respectively, for the two models and the GANIL condi- 
tions. At 6 = 10 fm we get an excitation probability of 3.4% and 3.9% in the 




(17) 
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first order calculation, decreasing with b~^. For impact parameters larger than 
^ fa 200 fm the excitation becomes adiabatic and the excitation probability 
drops off exponentially. The excitation probability is small compared to 1 and 
the non-conservation of unitarity in the first order approximation will not affect 
the calculation of the cross section. 

The apparent reduction of the B(El)-value (due to higher order effects) is 
given by 

R=-^ (18) 

which is plotted in Fig. 2 as a function of the projectile velocity v for the two 
model wave functions described above. We assume a minimum impact parameter 
bmin — 10 fm corresponding to a grazing collision of the projectile and target. 
This will give an estimate of the largest possible effect to be expected from the 
higher order contributions. We obtain a reduction, depending on the particular 
model chosen, of 5.5% or 10.1% for the GANIL energy and 3.7% or 6.9% for the 
RIKEN energy. The range of the reduction R in the two models gives a feeling 
of the reliabiliy of the results. The use of more realistic wavefunctions is outside 
the scope of the present work. The reduction of the excitation probability for the 
I state by higher order effects will be accompanied by an increase of the cross 
section in the breakup channel. 

In the dipole approximation, the lowest order correction was of 3'^'^ order. 
E1-E2 excitation contributes already in second order. Nevertheless, it is smaller 



than the third order El-El-El correction. For low ^ we can estimate the ratio of 
the two amphtudes as 



E1-E2 amphtude 1 v Z^jj 



(2) 

0.13 (19) 



El-El-El amphtude Za c {Z^]\y 
for Z = 82, Z^ff = 4/121, and Z^]\ = 4/11 at the GANIL energy. Thus the 
dipole approximation is reasonable, at least for a first exploration. 

Possible higher order effects in the ^^Mg intermediate energy Coulomb Ex- 
citation [Q] can also be estimated. Assuming a rigid rotor model, high energy 



Coulomb excitation was calculated in the sudden approximation in Ref. 
The characteristic strength parameter is 

C-^9l (20) 



Using the B(E2)-value found in Ref. p we have (see eq. 6 of Ref. ||T0|], where a 
factor TT is missing on the rhs.) 

1 fSvr 

Ql = —B{E2,0+^2+) (21) 

and 

C ^ 0.37 (22) 

with 6 = 15 fm for an energy of 49.2 ■ A MeV. From Fig. la of Ref. |Ty] it can be 
seen that higher order effects are negligible for the value of the strength parameter 
C. This is in agreement with the result found in the coupled channel calculation 
of Ref. i. 
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In conclusion, we provide a framework to apply corrections to I'** order elec- 
tromagnetic excitation. It is appropriate for fast collisions. We constructed oper- 
ators which take the influence of intermediate states into account. This can lead 
to a great simplification as compared, e.g., to the coupled channels approach. In 
this approach, a set of states, considered to be relevant, has to be chosen with 
known electromagnetic matrix-elements. In the present approach, of course, the 
model dependence cannot be altogether avoided; it enters when the corresponding 
matrix-elements of the operator has to be calculated, or at least estimated. 

As an example we studied the excitation of the | state in ^^Be. The impor- 
tance of higher order effects for this case of an extremely loosely bound nucleus 
was established. The estimate for the reduction of the cross section can only 
partly explain the observed reduction of the B(El)-value in the GANIL exper- 
iment. However, for a final analysis, more accurate calculations with improved 
wave functions, including E2 and nuclear effects in the excitation, should be 
performed. 

We wish to thank P. G. Hansen and M. G. Saint-Laurent for stimulating 
discussions. 

Note added in revision: In the meantime we got to know about a coupled 
channel study of the ^^Be Coulomb excitation by C. A. Bertulani, L. F. Canto, 
and M.S. Hussein. They get very similar conclusions as compared to our findings. 
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Figure Captions 

Figure 1: Reduction factor r for the excitation probability as a function of the 
characteristic parameter z defined in Eq. |l^ for the two model wave functions (I: 
— 11:---). 

Figure 2: Reduction R of the total cross section for the excitation of the \ 
first excited state from the )^ ground state of ^^Be as a function of the projectile 
velocity v for the two model wave functions (I: , II: ). 
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